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ABSTRACT 

Based on recent results on the frequency of Mgii absorption line systems in the "QSO 
behind RCS clusters" survey (QbC), we analyse the effects of the cluster environment 
on the sizes of baryonic haloes around galaxies. We use two independent models, i) 
an empirical halo occupation model which fits current measurements of the clustering 
and luminosity function of galaxies at low and high redshifts, and ii) the GALFORM 
semi-analytic model of galaxy formation, which follows the evolution of the galaxy 
population from first principles, adjusted to match the statistics of low and high 
redshift galaxies. In both models we constrain the Mgii halo sizes of field and cluster 
galaxies using observational results on the observed Mgii statistics. Our results for 
the field are in good agreement with previous works, indicating a typical Mgii halo 
size of TMgii — 50hy]^kpc in the semi-analytic model, and slightly lower in the halo 
occupation number approach. For the cluster environment, we find that both models 
require a median Mgii halo size of VMgii < lOhyj^'^kpc in order to reproduce the 
observed statistics on absorption line systems in clusters of galaxies. Based on the 
Chen & Tinker (2008) result that stronger systems occur closer to the Mgii halo 
centre, we find that strong absorption systems in clusters of galaxies occur at roughly 
a fixed fraction of the cold-warm halo size out to Ih^j^'^Mpc from the cluster centres. In 
contrast, weaker absorption systems appear to occur at progressively shorter relative 
fractions of this halo as the distance to the cluster centre decreases. These results 
reinforce our conclusions from Paper I and provide additional independent support 
for the stripping scenario of the cold gas of galaxies in massive clusters by the hot 
intracluster gas, e.g., as seen from X-ray data. 
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1 INTRODUCTION 

Understanding the influence of environment on galaxy evolu- 
tion either in the field or in galaxy clusters, provides impor- 
tant missing clues on the interplay of internal and external 
mechanisms that shape the galaxy population. In this work 
we concentrate on the statistics of Mg ll absorption systems 
in galaxies associated with high-redshift clusters, as recently 
obtained by Lopez et al. (2008; hereafter 'Paper I'). Our aim 
is to confront those statistics with different galaxy models 



in order to constrain the sizes of the cold-warm baryonic 
component of cluster galaxie^. 

Ever since the first studies of QSO absorption lines, 
Mgll was recognized as an excellent tracer of high-redshift 
galaxies (Bergeron & Stasinska 1986; Petitjean & Bergeron 
1990; Steidel & Sargent 2002; Churchill et al. 2000; Lanzetta, 
Turnshek & Wolfe, 1987, Tytler et al., 1987, Steidel & Sar- 
gent, 1992). These seminal works were then followed by nu- 
merous studies (Churchill et al. 1999; Nestor, Turnshek & 
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^ The temperature of the Mgn gas is expected to be a few lO** 
K, and therefore is termed here cold-warm gas. Notice that in 
some other cases, e.g., in Chen & Tinker (2008), Mgn is referred 
to as cold gas. 
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Rao. 2005; Nestor, Turnshek & Rao 2006, Narayanan et al. 
2007, Lynch, Charlton & Kim 2006; Prochter, Prochaska, 
& Buries 2006) that were able to distinguish between weak 
and strong Mgll absorption systems as possibly two dis- 
tinct populations. When selecting the stronger Mgll systems 
(with rest-frame equivalent width Wq'^'^^ > 0.3 A) the asso- 
ciated galaxies present a slight preference for blue, starburst 
galaxies with high metallicities (Zibetti et al., 2007, Elli- 
son, Kewley & Mallen-Ornelas, 2005). However, it is still not 
clear whether the stronger systems really occur in extended 
(hundreds of kpc) and virialised halos already in place at 
z = 2 (Churchill et al. 2005; Steidel et al. 2002) or if some 
are the signature of violent galaxy-scale outflows (Bouche et 
al. 2006; Prochter et al. 2006). 

Despite these extensive absorber-galaxy studies, and 
perhaps due to our lack of knowledge about the nature of 
Mg II absorbers, it has been difficult so far to establish corre- 
lations between absorption-line strength and absober-galaxy 
luminosity and impact parameter. Although some tentative 
correlations were found in the past (Churchill et al., 1999; 
Steidel et al. 1995), it is now becoming clear from deep imag- 
ing of the QSO fields, that not only the Mgll covering frac- 
tion must be less than unity (Chen & Tinker, 2008, Kacprzak 
et al., 2008), contrary to what was assumed for years, but 
also that there must be correlations with the galaxy envi- 
ronment. 

In Paper I we proposed to study the environmental de- 
pendence of the size of the cold-warm baryonic component 
using the frequency of Mgll in the line-of-sight (LOS) to 
QSOs close in projection to foreground clusters. Thus, con- 
trary to previous studies, our samples were constructed in 
an unbiased way by selecting a set of known galaxy over- 
densities and then searching for absorption systems in their 
neighbourhoods. In these particular lines of sight we ob- 
tained a comparatively fiatter equivalent width distribution 
with respect to what has been obtained for the field. This re- 
sult suggests that the cluster environment could be pruning 
the baryonic halo size significantly. 

Generally, the interpretation of the statistics on Mgll 
absorption systems in the fieltfl has been done following an- 
alytic approaches (as in Churchill et al., 1999). Only more 
recent works (Tinker & Chen, 2008, Chen & Tinker, 2008) 
use a model motivated by observational data within a cosmo- 
logical framework, the Halo Occupation Model (HO, Seljak, 
2000, Scoccimarro et al. 2001; Berlind & Weinberg 2002; see 
Zheng & Weinberg 2007, and references therein), to study 
the origin of the Mgll absorption systems, their geometry, 
and the dependence of equivalent width with the distance 
from a Mgll halo centre Q 



^ In this paper we refer to "the field" as a collection of objects 
which has been drawn at random from the full galaxy population, 
and therefore includes many different environments, from very un- 
derdense regions to high mass clusters. This coincides with the 
selection of QSOs made by several authors to detect Mg n absorp- 
tion systems, and will be adopted in our analysis of the numerical 
models. 

^ The HO models which we refer to in this paper, are rather com- 
plex Monte-Carlo simulations that reproduce density profiles of 
haloes of different sizes via fits to numerical simulation results, 
and populate these haloes with baryonic components following 
observational measurements of statistics such as the luminosity 



Building upon the results presented in Paper I, we per- 
form a theoretical study of Mg ll absorption systems using 
models to determine how the size of the cold- warm baryonic 
halo depends on environment. In order to do this we use 
two different models, (i) a HO model, constructed empiri- 
cally by Cooray (2006), adopting a ACDM cosmology, and 
(ii) a semi-analytic model by Bower et al. (2006, a version 
of the GALFORM model by Cole et al., 2000) which follows 
the evolution of galaxies in a ACDM Universe from first 
principles, tuned to reproduce the observed galaxy popula- 
tion at a wide range of redshifts (other semi-analytic models 
include Baugh et al., 2005, Croton et al., 2006, Lagos et al., 
2008, among others). We then concentrate on determining 
the cold-warm baryonic halo sizes and their dependence on 
galaxy environment, without modeling the distribution of 
gas in a direct way, but instead establishing the required 
cold- warm halo size needed to reproduce the observed Mg ll 
absorber statistics presented in Paper I for clusters of galax- 
ies. For simplicity, we start assuming a unity covering frac- 
tion and then discuss the implications of this choice. 

Throughout this paper the cosmological model adopted 
is characterised by the "concordance" parameters in line 
with estimates from the large-scale distribution of galax- 
ies and the temperature fluctuations in the cosmic mi- 
crowave background (Sanchez et al., 2006), namely, f2m = 
0.25, = 0.75, = 0.045, a Hubble constant Ho = 
71xh7ikm/s/Mpc, with hyi = 1, n = 1 and erg = 0.9. 

The paper is organised as follows. In Section 2 we briefly 
review the observational results from Paper I and the addi- 
tional reflnements applied to them to take into account the 
contamination from the fleld and the large-scale structure. 
Section 3 presents the two models in detail. Section 4 de- 
scribes the results, and a discussion is presented in Section 
5. Our conclusions are summarised in Section 6. 



2 OBSERVATIONS 

In this section we describe the available datasets used for the 
present work and the corrections we have made to account 
for contamination from the fleld and large-scale structures. 

2.1 The data 

In Paper I we used the Red sequence Cluster Survey (RCS, 
Gladders & Yee, 2005) to perform a search of clusters in 
the LOS to background quasars, and constructed two differ- 
ent sets of cluster-QSO pairs. The first set includes 46 high- 
resolution QSO spectra, and the second comprises lower res- 
olution spectra (these numbers do not include restrictions 
on signal-to- noise ratios). In all cases, the maximum physi- 
cal distance between the line-of-sight to QSOs and the clus- 
ter centres was set to 2h7^^Mpc, for a total sample of 529 
cluster-QSO pairs. 

We estimate the median mass of RCS clusters in the 
sample using the BgcR parameter (Gilbank et al., 2007) 
which can be used to estimate individual cluster masses. 

and correlation functions. All other non-HO models assuming sta- 
tistical distributions for observational quantities will be referred 
to simply as Monte-Carlo simulations. 
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For this work we only take into account results from clus- 
ters with at least a 2 — a detection of BgcR (that is, 2 — a 
away from zero). For this subsample of RCS clusters (all of 
which are paired to a background QSO for a total of 212 
pairs) we find a median mass of (l-64to 345) x lO^'^hfj^^Mo 
(errors correspond to the 20 and 80 percentiles of the mass 
distribution) which we will use in the remainder of this work. 

Clusters in the RCS survey are identified using a likeli- 
hood method which fits simultaneously for a projected clus- 
ter density profile and a red sequence at different rcdshifts 
(see Paper I for more details); this method produces an es- 
timate of the cluster photometric redshift as a by-product. 
A Mg II system is associated to a RCS cluster when the red- 
shift of the absorption falls within a 1 — cr interval around 
the cluster photometric redshift. In some cases a given ab- 
sorption system can satisfy this criterion for more than one 
cluster, in which case we test two different options, to assign 
the absorption system (i) to the closest cluster in projection 
(referred to as "nearest cluster", and (ii) to the farthest clus- 
ter in projection ("farthest cluster"). We will show that this 
choice does not affect our results in a significant way. respec- 
tively. 

Our subsamples arc defined via lower limits of equiva- 
lent width (EW). This is particularly important since with- 
out a clear modeling of the dependence of the line strength in 
our models we cannot exclude the occurrence of strong lines, 
particularly when taking into account indications that these 
occur closer to the centre of the absorbing galaxy (Chen 
& Tinker, 2008). We construct (i) sample S:WS (weak and 
strong systems) which includes all QSO-Cluster pairs from 
the high resolution spectra, and (ii) sample S:St (strong sys- 
tems only) which consists of QSO-Cluster pairs selected from 
the low- resolution data. Notice that sample S:St includes the 
QSO-Cluster pairs from S:WS (all the available pairs), but 
only considers absorption systems with EW above the lower 
limit of lA. We define hits in the same fashion as in Paper 
I, i.e., absorption systems within ±Az of a cluster redshift, 
where Az is the uncertainty in the cluster redshift. 

Notice that in order to maximize the number of hits, we 
do not impose the constrains on EW completeness or on pair 
redshifts used in paper I to correct for cluster completeness. 
As we show below, these more relaxed sample definitions 
introduce variations in our results that arc negligible when 
compared to the field and clustering corrections. 



2.2 



Correcting for systematic effects in the 
Cluster-QSO pair information 



The number of hits, Nhits, varies with the impact parameter 
to the cluster centre (see Tables 1 and 2). The number of 
QSO-cluster pairs is larger than the number of QSOs with 
measured spectra since in many cases there arc two or more 
clusters in the QSO line-of-sight. Given that a Mgll absorp- 
tion is associated to a cluster only if the redshift difference 
is within the cluster photometric redshift error, Az, there 
is a chance that some of these hits will have actually taken 
place in the field (which corresponds to the average environ- 
ment including voids and clusters, statistically) rather than 
in the associated RCS cluster. Therefore, we make a first 
correction to the number of hits by setting, 



where the correction from the field takes the form, 

NjiM = ^(«) X 2A2, (2) 

with z being the photometric redshift of the cluster, 2Az is 
summed over all clusters considered in a given sample, and 
dN/dz is the field estimate from the literature. For S:WS, 
dN/dz is given by 



^ = 1.9306 + iV,(l + z)" exp (-^(1 + «)"") , 

where the first term corresponds to the Churchill et al. 
(1999) estimate extended to 0.015A< VKq^™^ < O.SA, and 
the second term corresponds to the Nestor, Turnshek & Rao 
(2005) estimate for absorbers with EWs higher than Wo, 
which in this case is set to Wq = O.SA; their published pa- 
rameter values are iV, = 1.001(±0.132), a = 0.226(±0.170), 
W, = 0.443(±0.032), f3 = 0.634(±0.097) (uncertainties are 
shown only to illustrate the accuracy of the fit). For sample 
S:St dN/dz is given by the second term with Wo = lA. The 
adopted redshift corresponds to the median redshift of the 
RCS cluster sample, median(2:) = 0.6. The field corrections 
are listed in column 3 of Tables 1 and 2. 

It is well known that clusters occupy biased density 
peaks in the distribution of matter which are characterised 
by a high amplitude cluster-mass cross-correlation function 
(see for instance. Croft et al., 1997, Padilla et al., 2001). 
Hence, a correction has to be introduced to take into ac- 
count the enhanced matter density axound clusters which 
would increase the occurrence of Mg ii absorbers from the 
field in the surrounding few Mpc around each cluster. 

Wo then proceed to calculate the effect of clustering in 
the region surrounding the clusters, for which we propose a 
correction of the form, 



(3) 



Nhits 



N 
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■N^c 



(4) 



where N^corr includes the expected excess of field hits due 
to the ovcrdensity around our cluster sample. 



N^corr = I X Nfield 



(5) 



where the factor / is calculated using the cluster-mass cross- 
correlation function as 



/ 
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n {b,y'] 
(n) 



dy 



dy+ dy 



-1, 



(6) 



^hits = ^hits - Nfield, 



(1) 



where t he como ving distance in the redshift path is 
Yc = ''90 ~ bi^ ) '"cm indicates comoving distance, rgo = 
1.92h7j^Mpc is the median of the radius containing 90 per- 
cent of the cluster galaxies for our sample of RCS clusters, 
(n) is the average galaxy density, and n{b, y') is the galaxy 
density calculated using 

n{b,y') = {n)[l + ^Mib,y')]. (7) 

In the last expression, ^M{b,y') is the cluster-mass cross- 
correlation function, 

U{x,y') = ^MiVx^ + y'2) = $M(r) =P{^)\ (8) 

where ro = 5h~^Mpc and 7 = —1.8 parametrise the cor- 
relation function of the mass (see for instance Padilla et 
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range in Nf^n^ Field corrected Clustering corr. Number of QSO-Cluster pairs 

0.0 to 0.5 8 4.79 4.22 11 

0.5 to 1.0 5 2.41 2.05 12 

1.0 to 1.5 1 0.47 0.37 13 

1.5 to 2.0 1 0.48 0.40 10 

Table 1. Number of Mgll doublet hits in the S:WS sample (second column) for different ranges of impact parameter to the cluster 
centre (indicated in the first column). The third column contains the number of hits corrected from the contribution from the Field; the 
fourth column includes the effects of clustering for the correction; the fifth column shows the total number of QSO-Cluster pairs. 



range in 6/h^-^^Mpc AT^its Field corrected Clustering corr. Number of QSO-Cluster pairs 

0.0 to 0.5 1 0.96 0.95 14 

0.5 to 1.0 7 6.73 6.70 32 

1.0 to 1.5 2 1.92 1.91 55 

1.5 to 2.0 1 0.95 0.94 66 

Table 2. Number of Mgii doublet hits in the S:St sample. Columns are as in Table 1. 



al., 2004), and (3 = 1.6 represents the bias factor of haloes 
with the same median mass as our sample of RCS clusters 
(notice the change of notation for the bias with respect to 
previous works), obtained using the "concordance cosmol- 
ogy" and the Sheth, Mo & Tormen (2001) formalism. As we 
use the cross-correlation function between clusters and the 
meiss, there is only one factor of the bias parameter instead 
of the square of the bias as it would be the case in the cluster 
auto-correlation function. 

The eflect of these two corrections can be appreciated 
in the number of absorption systems associated to clusters, 
for different ranges of impact parameter with respect to the 
cluster centre shown in Tables 1 and 2. Notice that the effect 
of the corrections is very small for S:St, reinforcing the re- 
sults from Paper I where the frequency of strong absorption 
systems in clusters of galaxies was found to be significantly 
larger than in the field. 



3.1 Notation and definitions 

In the models we will make reference to two distinct types of 
objects, (i) cluster-mass dark-matter haloes, (ii) and galax- 
ies residing within these clusters. In all cases, the impact 
parameter will refer to the projected distance in physical 
units at the comoving distance to the cluster between the 
QSO line-of-sight and the centre of the cluster (real or simu- 
lated); we will use the notation b for this impact parameter. 
In the case of individual galaxies, we will refer to two sepa- 
rate quantities, (i) the scale- length of the gaseous disk, r^isk, 
and (ii) the extent of the Mg ll halo capable of producing a 
Mg II absorption feature above a minimum equivalent width, 
Wq''^^, that will be accordingly specified. The Mgii halo in 
the model galaxies is assumed to bo a uniform cloud of con- 
stant density with no holes (filling factor of unity; see Section 

5.2 for a discussion of the effects of a different value). Finally, 
in all cases the observational sample to be compared with 
the models will be characterised by its median equivalent 
width. 



3 MODELS 

We use theoretical predictions from numerical simulations to 

estimate the median size of the cold- warm gas halo from the 
available statistics on Mg ll absorption systems from cluster 
galaxies in the line-of-sight of QSOs. 

In order to do this we use two different methods. One 
consisting on a Monte-Carlo simulation and another on the 
analysis of the output of a semi-analytic model by Bower 
et al. (2006). The reason behind this choice is to allow our 
results to include a variance from the intrinsic assumptions 
present in very different modeling procedures, to allow a 
quantification of possible model-dependent effects. As men- 
tioned above, both models adopt the "concordance" cos- 
mology. In the following subsections we first present general 
definitions and notation, and then show the details of each 
of the models. 



3.2 Halo Occupation Monte-Carlo simulations 

The HO model is used in the literature to test whether the 
properties of individual galaxies depend on their host halo 
mass (e.g. Cooray, 2005); its parameters can be adjusted us- 
ing observed galaxy properties such as their clustering (sec 
for instance Zehavi et al., 2004) and luminosity functions in 
different bands (Cooray, 2006). A HO model does not in- 
clude physical prescriptions for galaxy evolution, although 
it can be used to infer the evolution of the population of 
galaxies in haloes via comparison with observations at vari- 
ous redshifts. 

In this first approach we use a HO model with pa- 
rameters adjusted by Cooray (2006) to match the lumi- 
nosity function of Sloan Digital Sky Survey (Abazajian et 
al., 2005) galaxies brighter than r-band absolute magni- 
tudes Mr = —17. This HO model assumes that the num- 
ber of galaxies per halo depends on the halo mass such 
that haloes with masses above a minimum value Mmin = 
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Figure 1. HO model: Number of events relative to the number of QSO-cluster pairs of single, double and triple hits (solid, dashed and 
dotted lines, respectively) as a function of median impact parameter measured with respect to the cluster centre. Left and right panels 
show the results from considering Vg^i = 50 and SOh^Tj^kpc, respectively. The cluster mass increases from the top to the bottom panels 
as indicated in the key. 



2.8 X lO^°h7/M0 host one central galaxy, and the num- 
ber of satellite galaxies per halo increases according to a 
power law Nsat = (M/Msat)^ , where M is the dark-matter 
mass of the host halo, Msat = IWtfl ^ lO^^hfj^A/o, and 
/? = 0.76 ±0.18. 

We combine this recipe with a projected Navarro, Frenk 
& White (1997, NFW) profile (from Yang et al. 2003) which 
we assume is followed by galaxies populating dark-matter 
haloes. The NFW profile allows us to determine the pro- 
jected density of galaxies as a function of the distance to 
the cluster centre, which is scaled according to the expected 
number of galaxies for clusters of a given mass, obtained 
from the HO model. The projection of the NFW profile is 
done by integrating the 3-dimensional profile on the direc- 
tion of the line-of-sight out to 2 virial radii from the cluster 
centre (see Yang et al., 2003, for more details). We use this 
projected galaxy number density, (j{b), as a function of the 
projected distance to the cluster centre to determine the 
median Mgll halo size, rMgii, necessary to produce a given 
frequency of hits, Nhits/Nhos, in a sample of Nlos QSO- 
cluster pairs using. 



where we have assumed that all galaxies have equal Mgll 
halo sizes. 

Figure [T] shows the expected number of hits per QSO- 
cluster pairs for clusters with median masses of lO^^'^h^^ Mq 
(top panels), increasing to lO^^'^h^^ Mq (bottom panels). 
We show the results for 1, 2 and 3 simultaneous hits (solid, 
dashed and dotted lines) in the LOS of one target quasar 
as a function of the median impact parameter, bm, calcu- 
lated at the median redshift of the RCS cluster sample, 
median(z) = 0.6. In this case, we assume two values for 
the median Mgll halo sizes, rMgii = 50 and SOh^j^kpc (left 
and right panels, respectively), which correspond to a range 
of sizes consistent with estimates from absorption systems 
of field galaxies in QSO spectra (Churchill, 2001, Churchill 
& Vogt, 2001, Nestor, Turnshek & Rao, 2005). The models 
show that it is easier to obtain more events of multiple hits 
when observing more massive clusters as well as when choos- 
ing background QSOs closer to the cluster centres. Also, a 
larger Mg ii halo will produce more hits, as expected. These 
plots allow us to study the particular case of LOS number 
14 from Paper I, where two Mgll absorption systems over- 
lap with the allowed ranges of photometric redshifts of 6 
individual RCS clusters. Using the expected frequency for 
single and double hits for the median mass of our RCS clus- 
ter sample, the resulting likelihood that both Mg ll systems 
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Figure 2. HO Model: Relative errors in r'Mgii as a function of 
number of target QSOs for a fixed cluster mass of M = 1.64 x 
IO^^Ii^^'^Mq and median impact parameter bm = ISSQh^j'^kpc. 
Results are shown for rjyigii = 56, 63 and 77h.^-^^kpc. The grey 
vertical line shows the number of QSO-cluster pairs in the S:St 
sample. 

correspond to only one out of the six clusters is ~ 10"' times 
lower than the case where the individual absorptions are as- 
sociated to two different clusters (out of the available 6). 
This calculation corresponds to the average impact parame- 
ter hm = 1200h^j^kpc, characterising this LOS, and consid- 
ers the measured median masses of the RCS clusters in our 
sample (in the case where the double hit and single hit are 
equally likely, the former are only ~ 5 times less likely to 
happen than the latter for this LOS). Therefore, from this 
point on we consider this case as two single hits (Table 1 is 
constructed accordingly). This would not be the case if the 
multiple absorption system was consistent with the photo- 
metric redsliift of a single RCS cluster. This particular LOS 
offers an excellent laboratory for follow-up spectroscopy of 
absorbers. 

We estimate the relative errors in Mgll halo sizes, 
^TMg II /rMg II, assuming Poisson statistics and using the 
frequency of single hits and samples of different total num- 
ber of QSO-cluster pairs. For this calculation we adopt fixed 
values of cluster mass, M = 1.64 x lO^^'h^^^Mo, and impact 
parameter, b = ISSOh^^^kpc, corresponding to the median 
values in our full sample of QSO-RCS cluster pairs. In Figure 
2 we show the variation of Arjwg ii /I'Mg ii for three different 
Mg II halo sizes as a function of the number of QSO-cluster 
pairs. The vertical grey line indicates the expected errors in 
the Mg II halo size for the total number of 529 pairs in sam- 
ple S:St, and indicates that for sizes larger than ~ 50h^^^kpc, 
the uncertainties would be below a 25 percent of the inferred 
size. 

We will use this model to infer the typical size of the 
Mgll halo and its associated measurement error in the ob- 
servational data described above. 

In what follows, we introduce the details of the addi- 
tional model which uses a cosmological numerical simula- 



tion in combination with a semi-analytic model of galaxy 
formation. 



3.3 GALFORM Semi-analytic galaxies 

This model differs from the HO in that it follows the evolu- 
tion of galaxies using a number of physical prescriptions for 
the different galaxy components, including cold gas, stars, 
metallicity of stars and gas (hot and cold), and the proper- 
ties of the central super-massive black hole. However, there 
is a strong observational component to the model since its 
parameters are tuned to match the observed galaxy popula- 
tion at various redshifts. 

The main advantage in using a semi-analytic model re- 
lics on the following. The global properties of large concen- 
trations of matter either in simulations or in observations 
are well known on average; however, the complex nature of 
the non-linear collapse that formed virialised objects has im- 
portant consequences on their diversity, even for objects of 
similar masses. The properties that are not included in a 
general description of clusters include their asphericity, sub- 
structure and formation history, and their relation to the 
surrounding environment. Therefore, the use of a cosmolog- 
ical numerical simulation allows us to include this diversity 
in the population of cluster-size dark-matter haloes, an effect 
that is very difficult to include in the model presented in the 
previous subsection, which relies on Monte-Carlo modeling. 

In this work we use the Bower et al. (2006) GAL- 
FORM galaxy catalogue which conforms a complete sample 
of galaxies down to a magnitude = —17. These galaxies 
populate the Millennium Simulation (Springel et al., 2005), 
which follows the evolution of ~ 1 x 10^° collisionless dark- 
matter particles from z = 50 to the present on a box of 
704.2hfi^Mpc a side, for a mass resolution per particle of 
1.21 X lO^hyj^M©. This simulation allows to identify bound 
dark-matter structures and to follow their descendants at 
different redshifts. This, in turn, can be used to infer the 
galaxy population within each dark-matter halo, via various 
assumptions regarding the different processes that shape the 
formation and evolution of galaxies in a semi-analytic model 
(Croton et al., 2007, Malbon et al., 2007, Lagos, Cora & 
Padilla, 2008, De Lucia et al., 2006). 

It is important to note that these galaxies do not include 
a description of their Mgll content, and therefore can only 
be used in a similar way to the HO model procedure. Given 
that the model follows the evolution of the disk scale-length, 
^disk of each individual galaxy with time, we assume that the 
radius of the Mg ii halos will be, 

rMgii = Ax Tdisk, (10) 

where A is a proportionality constant that will be varied in 
order to match the observed frequency of Mg ll absorption 
events. This is an important hypothesis that we will adopt 
throughout this work, namely, that the size of the Mg ii halo 
is proportional to the disk size. 

Our sample of simulated clusters is constructed by se- 
lecting all systems above a minimum mass such that the 
median mass equals that present in our RCS cluster sample. 
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3.3.1 Simulated composite cluster 

In order to study the frequency of galaxies in front of back- 
ground QSOs in the numerical model, we construct a com- 
posite cluster using all the objects within our sample of 
simulated clusters. We do this by stacking all the clusters 
using their most bound particle as their centres; the semi- 
analytic model by Bower et al. (2006) uses this position to 
place the central galaxy of each halo. The stacking proce- 
dure erases all cluster membership information. Therefore 
our simulated composite cluster cannot be used to study 
multiple absorption events from individual clusters. As was 
mentioned above, the multiple hit candidate present in the 
data from paper I is considered as two single hits due to its 
low probability of occurrence. 

Given that the population of background QSOs is at 
a very large distance behind the clusters in our sample, we 
assume that the former are distributed at random within 
2h^^^Mpc from the cluster centres, in projection. 
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3.3.2 Morphological types 

Among the information available for each galaxy in the semi- 
analytic model are the stellar mass in the disk and bulge 
components. Therefore, it is possible to separate a popula- 
tion of elliptical and spiral galaxies, which may prove im- 
portant in light of previous results suggesting the presence 
of strong Mg ll absorption systems mainly in spiral systems 
(Zibetti et al., 2007). 

Therefore, our analysis of number of galaxies in the line- 
of-sight of background QSOs is done to a first approximation 
only on spiral galaxies selected so that the stellar content in 
their bulges is only unto 70 percent of the total stellar mass 
(Bertone et al., 20070. 

We bear in mind that Zibetti et al. (2007) also show 
that there may be dependencies of the Mgll absorption on 
galaxy colour and morphology. Given the simplicity of the 
approach we follow in this work, we will not attempt to 
model this dependency but instead will offer estimates on 
the expected impact of the morphological selection on our 
results in Section [5.21 



4 RESULTS 

In this section we provide separately the results obtained us- 
ing the two models presented in the previous section, firstly 
from the observationally motivated Halo Occupation model, 
and secondly from the semi-analytic model galaxies in the 
Millennium simulation (Springel et al., 2005). 

4.1 HO Model 

We use Eq. |5] to determine the sizes of galaxy Mg ll haloes 
that would produce the rate of hits shown in Table 1 (only 
calculated for the S:WS sample), for clusters with median 

* Bertone et al. find that in the case of the semi-analytic model 
we are using in this work, a limit of 70 percent of the total mass 
in stars is the minimum fraction for early type objects that allows 
to reproduce the observed dependence of morphological fractions 
as a function of stellar mass and local density. 



Figure 3. HO Model: Mgll halo sizes as a function of impact pa- 
rameter, according to the HO model results. Open squares show 
the results for the total number of hits in the S:WS sample, open 
pentagons correspond to the number of hits corrected by the ex- 
pected field hits, and open triangles show the results from adding 
the correction for clustering. The open star symbols indicate the 
resulting sizes when one of the cluster members is positioned at 
the centre of mass of the cluster. The horizontal solid line indi- 
cates the radius of Mgn haloes in the field obtained from the HO 
model. 



mass 1.64 x IQ^'^h^l Mq. We show the resulting typical sizes 
in Figure |3] for the total number of hits in the S:WS sample 
(open squares) and also for the corrected counts by the ex- 
pected contamination from the field (open pentagons) and 
clustering (open triangles). The results when considering the 
"farthest cluster" variant of our S:WS sample fall within the 
errorbars of these estimates, and therefore are not shown. 
Errorbars are calculated by assuming Poisson statistics. The 
points along the abscissa corresponds to the middle of the 
bin in b, with small shifts for the different cases, added to 
improve clarity. 

An important change in the inferred Mg ll halo sizes is 
produced when one of the cluster members is placed at the 
centre of the cluster (i.e. a central galaxy). This affects the 
innermost region for which it is straight-forward to infer the 
expected fraction of line-of-sights that would pass through 
the central galaxy Mgn halo once its typical size, rMgii, is 
known. This is an iterative process which converges rapidly 
after three iterations. When taking this into account, the re- 
sults change quite dramatically for the low impact parameter 
bin as it is shown by the open star symbols in Figure[31 These 
results indicate that near the cluster centres, the Mg ll halo 
of a galaxy could either be as large as rMg ii — 16hY^^kpc, or 
as small as rugii — 2h^j^'^kpc, depending on whether a cen- 
tral galaxy is assumed to lie exactly at the observationally 
determined cluster centre (which is subject to uncertainties). 

Our results show that the typical size of the Mg ll halo 
tends to increase towards the outskirts of the cluster. This 
result is independent of the corrections applied to the obser- 
vational data, and is not sensitive to the presence of a central 
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Figure 4. Semi-analytic model: Histogram of galaxy disk scale- 
lengths (lower X-axis) and Mgll galaxy haloes (upper x-axis) 
needed to reproduce the number of hits reported for the field. 
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Figure 5. Semi-analytic model: Dependence on impact param- 
eter to the cluster centre of the median disk scale-length in the 
semi-analytic model. The errorbars show the error of the median. 



galaxy at the centre of mass of the cluster. This is also in rea- 
sonable agreement with previous studies that report that the 
cold/cold-warm gas haloes around galaxies may be stripped 
off by the hot intracluster gas (see for instance Heinz et al., 
2003). 

We also use the HO model to estimate the typical Mg ll 
halo size in the field, which is indicated by the horizon- 
tal grey solid line in Figure |3l We estimate this size by 
integrating over the full range of halo masses where dark- 
matter haloes are expected to have at least one galaxy with 
Mr < —17, for the EW range corresponding to sample 
S:WS. As can be seen, our estimate of r-Mgii — SOh^j'^kpc (in 



broad agreement with observational estimates, e.g. Churchill 
et al., 1999) is slightly larger than our findings for the inner 
cluster environment, in agreement with the stripping sce- 
nario in clusters. For reference, the galaxy number density 
in the HO model is n = 7.639 x lO'^/hfj^Mpc^. Notice that 
this result corresponds to the HO model; the results from 
the analysis of the semi-analytic model for the field are pre- 
sented in the following subsection. 

4.2 GALFORM semi-analytic model 

Given the better description of individual clusters allowed 
by the combination of a numerical simulation and a semi- 
analytic model, the results from the composite cluster may 
provide further clues regarding the dependence of Mg ll halo 
on the distance to the cluster centre. 

We start by studying the expected typical halo size in 
the field. In order to do this, we randomly select positions in 
the simulation and take all the galaxies within 2h^j^Mpc in 
projection from these centres stacked together. This stack is 
characterised by a comoving length which we convert to a 
redshift path. We then place one QSO in the background and 
sort all the galaxies in the stack with respect to their pro- 
jected distance to the QSO, and assume that their Mg ll halo 
sizes will be directly proportional to their disk scale-lengths, 
as was mentioned in Section 3.3. All the Mg ll haloes defined 
this way that contain the LOS to the QSO are counted. We 
then vary the constant A until the simulation reproduces 
the same counts per unit redshift path, dN/dz, reported by 
the study of the field by Churchill et al. (1999) 

Figure [3] shows the results of the semi-analytic ap- 
proach. It depicts the distribution of disk scale-lengths in the 
field (lower x-axis label), where galaxies show a wide range 
of values with a peak at rdisk = Gh^^^kpc. The top x-axis 
label shows the resulting Mgll halo sizes, with a distribu- 
tion peak at a value of rMgii — SOh^^^kpc, well in agreement 
with the observational estimate of Churchill et al. (1999) 
for EW> O.O2A. Also note that our estimate from this ap- 
proach is only slightly higher than the results from the HO 
presented in the previous section. This difference arises from 
the different EW limits assumed in each case (otherwise both 
results agree with each other). 

The analysis of the composite cluster constructed from 
the GALFORM model also provides useful insights on the 
size of the Mgll haloes in the cluster environment, and par- 
ticularly as a function of distance to the cluster centre. In the 
semi-analytic model, galaxy disk scale-lengths show smaller 
typical sizes near the cluster centres as can be seen in Figure 
[S] where we show the median size of galaxy disks (points; 
errorbars indicate the error of the mean assuming Poisson 
statistics). This is a first indication that the model indeed 
includes important environmental effects on the galaxy pop- 
ulation (notice that these values are not to be compared to 
the field Mgll region sizes denoted as rMgii). In the GAL- 
FORM semi-analytic model, the physical process driving the 
environmental changes within clusters is the removal of the 
hot gas reservoir around galaxies once they enter a new dark- 
matter halo; this process effectively stops the cooling of gas 
onto the galactic disk halting its growth. 

Following the procedure used to infer the typical sizes 
of Mg II haloes in the field, we also find the proportionality 
constant that multiplied by the disk scale-lengths gives a 
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Figure 6. Semi-analytic model: Histogram of galaxy disk scale-lengths (lower x-axes) and Mgll galaxy haloes (upper x-axes) needed 
to reproduce the number of hits in the observations, for different ranges of impact parameter b, shown in the key (increasing values from 
left to right panels) for the S:WS and S:St samples (top and bottom panels, respectively). 
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Figure 7. Semi-analytic model: Median sizes of Mg II haloes as 
a function of projected distance to the cluster centre. Errorbars 
show the error of the median. 



Mgll cloud size that provides a match to the abundance of 
Mgll absorption systems in cluster-QSO pairs in the S:WS 
and S:St samples. Figure [6] shows histograms of galaxy disk 
scale-lengths (lower x-axes) and resulting typical Mg ll halo 
sizes (top x-axes) for two ranges of impact parameter (in- 
creasing values of b from the left to the right panels) for 
the S:WS (top) and S:St (bottom) samples. The size dis- 
tributions show a marked raise towards a single peak (on 
the position of these peaks we will centre our discussion) 
and a slower fall-off for higher values of r^isk, with a dis- 
tribution width that clearly increases towards the cluster 
outskirts, in both disk and Mgll halo sizes. From the loca- 
tion of the peaks it can be seen that the increase of the disk 
scale-lengths with the impact parameter shown in Figure 5 
is confirmed. It can also be seen that the location of the 
Mgll size distribution peaks show a similar trend of larger 
values further away from the cluster centres. Furthermore, 
the Mgll halo sizes are consistent with larger values for the 
S:WS sample for the lowest impact parameter bin, which re- 
flects the fact that the frequency of Mg ll absorption events 
is higher in this sample than in S:St at such distances from 
the cluster centre. 

Figure [7] shows the median size of a Mgll cloud as 
a function of the distance to the cluster centre (errorbars 
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show the error of the median) for the S:WS (triangles) and 
S:St (circles) samples. As can be seen from both figures, the 
trend of a larger galaxy size at larger distances from the 
cluster centre is clearly found in both samples with values 
ranging from rMgii — 2 — Gh^^^kpc at 6 < O.Sh^j^'^Mpc to 
rMgii — 25h7^^kpc for b > l.Shy^^Mpc, in good agreement 
with the HO results with a central galaxy[f|. The correspond- 
ing inferred size for typical Mgll haloes in the field for the 
range of EW in sample S:WS for the semi-analytic model 
is 40h^^^kpc (notice that this value differs slightly with the 
result from Figure 4, due to the different EWs considered). 

We remind the reader that these results are motivated 
by the observational frequency of Mg ll absorption systems 
in our samples of RCS cluster-QSO pairs (the proportion- 
ality constant in Eq. 10 depends on the impact parameter). 
Therefore, the trends of smaller Mg ll halo (cf. Fig. (Tjl and 
disk sizes (cf. Fig. [Sjl towards the cluster centres are inde- 
pendent results. Furthermore, the reasonable consistency be- 
tween the results from the HO model and the semi-analytic 
simulation provides further reliability to the findings on the 
trends of Mgll halo size with the distance to the cluster 
centre. 



5 DISCUSSION 

In the previous section we studied the inferred typical Mg ll 
halo sizes using two different models, the HO and the output 
from a semi-analytic simulation. We find that both mod- 
els are able to reproduce the observed frequency of Mgll 
absorption events found in observational data in Paper I. 
More importantly, our results indicate that the Mgll halo 
sizes tend to increase from relatively low values of rjvig ii — 
2 - 6hf;^^kpc to rMgii ~ 26 - 40hf;^^kpc for b = O.Shf^^^Mpc 
to 6 > l.Sh^j^'^Mpc; both models produce consistent results. 

In what follows, we present a study of the fraction of the 
Mgll halo responsible for producing different Mgll absorp- 
tion line strengths, and an analysis of possible systematic 
effects present in our current measurements. 



5.1 Regions of strong Mgll absorption line 
systems 

We take full advantage of the additional information avail- 
able in the semi-analytic model regarding the galaxy disk 
scale- lengths, and combine this with previous results by 
Chen & Tinker (2008) who find that stronger absorption 
systems are produced by clouds closer to the galaxy cen- 
treS 

to study the typical size of the regions in the Mgll 
halo responsible for different absorption line strengths. 

In our analysis of the semi-analytic model we maintain 
a record of the distance, d, between each QSO line-of-sight 



^ Notice that the Mgn halo size at 6 > 1.5hy-^"'^Mpc is not 
expected to resemble the Mgn halo size inferred for the field 
since this sample contains only galaxies in high mass dark-matter 
haloes, and therefore does not include a representative sample of 
field galaxies, even in the cluster outskirts. 

^ In their terminology this would be a halo centre, which can be 
applied to each galaxy in the semi-analytic model since these arc 
located at the centres of bound substructures (or sub-haloes) in 
the numerical simulation. 
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Figure 8. Semi-analytic model: Extent of different absorption 
line equivalent widths in units of the Mgll halo radius, as a func- 
tion of impact parameter to the cluster centres. The results from 
the S:WS sample are shown by triangles, and from the S;St sam- 
ple by circles. Open symbols correspond to strong absorption line 
systems; filled symbols to weak absorption lines. 



and each galaxy centre. Therefore, after finding the propor- 
tionality constant A (Eq. I10|l that reproduces the frequency 
of hits found in Paper I, we can sort the galaxies with respect 
to d. Then, for different ranges of cluster-centric impact pa- 
rameter, we study the observed fractions of different line- 
strengths, and find the galacto-centric distance out to which 
the same fractions of model galaxies are included. Such a 
distance is then identified as the maximum distance out to 
which a given line-strength would be produced in the model. 
Figure [S] presents the resulting fraction of the Mg ll halo size 
out to which a given line-strength occurs, as a function of 
distance to the cluster centre. This rescaling is adopted in 
order to factor out the measured variation of the global Mg ll 
halo size with the impact parameter to the cluster centre. 
Open and solid symbols represent high [Wq"^^^ > 2A) and 
low (Wq'^^'' < 2A) equivalent width absorption systems, re- 
spectively; triangles correspond to results from the S:WS 
sample, circles to S:St. In the case of the high equivalent 
width results, the large errorbars only allow a possible re- 
jection of an increase of the Mgll halo region producing 
such absorption events. However, there is a mild indication 
that the low equivalent width absorption systems (all from 
S:WS) tend to occur at increasingly larger distances from 
the galaxy centre in terms of the Mg ll halo size towards the 
outskirts of clusters. These results are again in good agree- 
ment with the stripping scenario, where the low density Mg ll 
clouds within a given galaxy would be more easily removed 
than the denser Mg ll regions nearer the galaxy centre that 
would be responsible for high equivalent width absorption 
systems. 
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5.2 Possible systematic effects 

The parametrisation of both models, HO and semi-analytic, 
used to mimic the observational procedure carried out in 
Paper I, depends on a number of parameters, assumptions 
and measurements that can induce systematic effects. 

For instance, the results shown in the previous sections 
were obtained for clusters of a median mass corresponding 
to the estimates for the RCS clusters. In order to assess the 
possible systematic biases arising from errors in these clus- 
ter mass measurements we change the adopted median mass 
by factors of 4 and 0.25, approximately 0.6dex in the clus- 
ter mass. The results in the HO and semi-analytic models 
are affected by very similar relative variations, which do not 
show a dependence with the distance to the cluster centre. 
When the lower value for the median mass is adopted, the 
inferred median sizes of the Mgll haloes increases by ~ 10 
percent; alternatively, a higher median mass produces a de- 
crease in the inferred size by a 20 — 25 percent. Given that 
the effect does not vary with the distance to the cluster cen- 
tre, our conclusions would not be significantly affected by 
this source of systematic uncertainties. 

It should also be taken into account that in our model- 
ing, we have considered spherical Mgll haloes. Given the 
lack of detailed information on the thickness of cold gas 
disks, we do not venture into assuming a given disk thickness 
in this work, as has also not been done in previous works 
studying the size of Mgll absorption regions in the field as 
in ChurchiU et al. (2001). 

It has been argued that a unity covering fraction for 
Mgll haloes is rather unlikely (Bechtold & EUingson, 1992, 
Tripp & Bowen, 2005, Chen & Tinker, 2008, Kacprzak et 
al., 2008). Since the models we have adopted here assume a 
fixed value for the covering fraction (which we have chosen 
to be C = 1), we cannot study a possible dependence with 
environment. However, it is easy to see that a different cov- 
ering fraction will change our results on the Mgll halo sizes 
according to rMgii to r'iign = -^rMgij. 

Finally, the effect from considering absorption from low 
gas content galaxies should also be taken into account. In 
the semi-analytic model, the fraction of elliptical galaxies 
changes as a function of the distance to the cluster centre. 
This is specially true for r/r^ir < 1, where the fraction of 
model elliptical galaxies (see Section I3.3.2p increases from 
~ 15 percent to ~ 50 percent at rjr^ir = 0.015, measured 
in projection (these fractions are in good agreement with 
results from the SDSS from Goto et al., 2003). At larger 
distances from the cluster centres, the fraction of elliptical 
galaxies remains practically unchanged. For sample S:WS 
(and similarly for sample S:St), considering an absorption 
from elliptical galaxies in the semi-analytic model would in- 
crease by a factor of ~ 1.43 III) the number of potential ab- 
sorbers and therefore decrease our estimate of the typical 
Mgll size inside clusters by a 16.4 percent. Our estimates 
at larger distances from the cluster centre will be equally 
affected by a decrease of approximately 7.8 percent. The 

^ This is an upper limit for the contribution of ellipticals to the 
projected galaxy number density at a median impact parameter 
of b/r„ir = 0.17. This value corresponds to the first bin in our 
measurements of TMgii, b < SOOh^j^kpc, for the median virial 
radius of the clusters, r^i^ = 775h^j"'^kpc. 



latter estimate would also change our estimates for the typ- 
ical Mg II halo size inferred for the field accordingly in both, 
the semi-analytic and HO models. In consequence, our con- 
clusions do not change significantly due to this effect. 



6 SUMMARY AND CONCLUSIONS 

In this paper we have presented a theoretical study of Mg ll 
halo sizes of galaxies in clusters and in the field, by com- 
paring recent observational results on the incidence of Mg ll 
absorption systems in RCS cluster-QSO pairs by Lopez et 
al. (2008), and the predicted systems from the Halo Occu- 
pation and Semi-analytic models. 

As a first step in our analysis we applied corrections 
to the observed frequency of absorption events to take into 
account possible contaminations from the field, and effects 
of the Large-Scale Structure. 

The results from the HO and semi-analytic models indi- 
cate that the typical field Mgll halo is rgai ~ 39 — SOh^j^'^kpc 
(for EW> O.O2A), a value in broad agreement with previous 
results pointing towards rgai ~ 50 — lOOh^^^kpc found for 
similar absorber EWs (Churchill et al., 2001, Zibetti et al., 
2007, Chen & Tinker, 2008, Kacprzak et al., 2008). 

The HO model also indicates that the Mg 11 halo tends 
to be smaller near the cluster centres, reaching values of 
Tgai — 16h^^^kpc. When one of the cluster members occupies 
the centre of mass of the cluster, the Mg 11 halo size can be as 
low as ~ 2h7^^kpc. This result is obtained assuming a fixed 
galaxy size for all the cluster members. The results from the 
Semi-analytic model are consistent with those from the HO 
showing a clear decrement in the typical Mgll size towards 
the cluster centres. An interesting by-product of this analysis 
is that the galaxy disk scale-lengths in the model also show a 
clear increase towards the outskirts of clusters, in accordance 
with the scenario where strong interactions strip galaxies of 
their outer disks in cluster environments. 

We also used the semi-analytic model to estimate the 
radii of Mg 11 haloes producing different strengths of absorp- 
tion. In order to do this we assumed that stronger lines are 
produced nearer the centres of galaxies (following Chen & 
Tinker, 2008). We found that in terms of the Mgll halo 
size, the strong absorption regions occur out to a rather 
fixed value, independent of the distance to the cluster centre. 
However, weaker absorption line systems show a tendency 
to occur at larger distances from the galaxy centre (in units 
of their global Mgll halo radius) towards the outskirts of 
clusters. This trend could explain the fiatter EW distribu- 
tion found in paper I, that included all impact parameters 
within b < 2h7i^Mpc. 

Our results indicate that the effect of the cluster envi- 
ronment on galaxies is extremely important and can produce 
a decrement of up to a 90 percent of the Mg 11 halo size with 
respect to galaxies in more gentle environments such as the 
field. 
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